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ANALYTICAL SOLUTIONS OF SPECIFIC CLASSES OF ASTROPHYSICAL 

RADIATING SHOCKS 

E. Falize^'^, C. Michaut^, S. Bouquet^'^ and C. Cavet^ 

Abstract. In this paper we study specific classes of radiating sliocks wfiidi are widely spread in astro- 
physical environments. We present more general solutions of their structure and proceed to the analytical 
determination of physical quantities. 

1 Introduction 



Radiating shocks play a crucial role in astrophysical environments (Mignone 2005) as well as in laboratory 
plasmas generated by powerful facilities (Drake 2006). In astrophysics, we can find them in the head of stellar 
and galactic jets (Blondin & CiofR 1989), in the first phase of supernovae explosions, in the late phase of 
supernova remnants (Chevalier & Blondin 1995), in accreting systems such as the magnetospherical accretion 
of T Tauri (Gunther et al. 2007) or magnetic cataclysmic variables (Cropper 1990). These shocks are at the 
basis of several models which allow the interpretation of astronomical observations. For instance, they may 
explain the recent observations of X spectra (Gunther et al. 2007) of classical T Tauri: TW Hya, BP Tau and 
V4056Sgr. Thus, it is very important to understand the structure and stability of these shocks. The structure of 
radiating shocks depends on the nature of pre-shock and post-shock media (Drake 2006). In this paper, we will 
consider shocks for which radiative effects can be modelled by entropy losses. This is the more common shock 
that it is encountered in interstellar phenomena. We will derive a solution by assuming a power law model for 
the cooling function, which generalizes the five specific solutions we find in the literature (Chevalier & Imamura 
1982). Firstly, we will recall the theoretical model and discuss the approximation. Secondly, we will present 
the analytical solution for the one- and two-cooling processes problems. Finally, we proceed to the evaluation 
of fundamental physical quantities that characterize this type of shock and compare them with the results we 
already have. 



2 Approximations and theoretical post-shock medium modelling 

We consider a plane-parallel [d/dy = d/dz = 0) collisional shock (i.e., tu « tdyn where ta is the characteristic 
time of collision between ions and tdyn is the dynamical time) with a post-shock medium which can be defined 
by a single temperature model {i.e., tei << tcooi where td is the characteristic time of energy exchanged between 
electrons and ions and tcooi is the radiative cooling time). Moreover, we suppose that the local gravitational 
field does not modify the shock structure {i.e., g*Xs << v'^ where g*, Xs and Vs are respectively the gravitational 
field, the thickness of the cooling layer and the shock velocity). Finally, we assume that the shock is stationnary 
(i.e., tcooi << tdyn which means that cooling effects are faster than dynamical ones). Thus, the equations that 
give the evolution of post-shock medium are (Kylafis & Lamb 1982): 



^H=0 -^[pv^+P]=0 
ax ax 



dP 
dx 



Pdp 
^ p dx 



-{-f-l)A{p,P,x) 



(2.1) 



where x, p, v, P, 7, A are respectively the spatial coordinate, the density, the velocity, the pressure, the 
polytropic index and the cooling function. Although these equations are theoretically consistent, we must 
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specify a equation of state in order to connect the microscopic phenomena to the cooHng function. We consider: 
P = e[Z\p"'T^ where T, e[Z\, a and fi are respectivefy the temperature, a function of ionization Z and two free 
exponants where we must impose 7(1 — (3) = {a — 13) in order to make sure that the entropy is preserved in the 
post-shock region. On the shock front we have the Rankine-Hugoniot relations which are satisfied: 
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Fig. 1. Schematic representation of radiating shock with the definition of the different physical quantities and notations. 
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(2.2) 



where po and M are respectively the pre-shock density and the Mach number (see the figure 1). It appears that 
a natural rescaling is the normalisation of spatial coordinate to the thickness of the cooling zone (^ = x/xs), the 
density in the post-shock region to the pre-shock one, the velocity to the shock velocity and the pressure to the 
ram pressure. In order to reduce the number of equations (j2.ip we introduce an intermediate function ri[^] defined 
by v[S,] = —Vsr][£^]. The first and second equations of (|2.ip imply = Po/vi^] ^^^d P[^] = PoV'^{0'iv[S,] + 02) 
where oi and 02 are two constants. This transformation is consistent with the Rankine-Hugoniot relations 
provided ry[l] = [2 -I- (7 - + 1)M'^] and P[^] = po«?(l + l/^M^ - r][(]). The determination of the 

intermediate function is possible through the last equation of (j2.ip which is given by: 



(7 + l)ry-7 1 + 



7M2 



drj 
d^ 



^(7-i)-^A(p,p,e) 



Pov^ 



(2.3) 



It is this equation that we examine now. 



3 Analytical solutions to one- and two-power laws cooling functions problem 

First, we suppose that the cooling function writes: A[p,P,x] = Aqp'^P''{x + Xq)^ where Aq, e, (, 9 and Xo 
are five constants. This form generalizes the optically thin case where 9 — 0. The arbitrarity of e and C is 
motivated by the fact that A oc KpaT^ where Kp is the Planck opacity which can be modelled by a power law at 
high temperature. Introducing this form in equation (|2.3p and using the integral form of Gauss hypergeometric 
function (that we note F2i[a, 6; c; a;]), we obtain an implicit form of the compression profile: 



cste + (7 + 1)?7' 



e+2 



1 



-c 



77? 



ni-C 



7M2 



r[e + 2] 



F21 



C,e+l;e + 2 



r[£ + 2] 

r[e + 3] 



^21 



C,e + 2;e-H3 



1 + 1/[7M2] 



Koln[^ + Xo/xs]; 9 = -I 



(3.1) 
(3.2) 



Analytical Solutions of Astrophysical Radiating Shocks 



3 



where kq = (7 — l)x^'^^ Aopl'^''~^v'^''~^ and cste is defined by a specific value of the implicit function. This 
general solution is an extension of the five specific solutions known in the literature (Chevalier & Imamura 
1982). Now we consider the case with two optically thin cooling processes. The cooling function writes: 
A(/9, P) = Ao^ip'^iP^i + Aq^2P'^^P''^ where we suppose, without loss in generahty, that tcooi.i < tcooi.2- Thus, the 
development of inverse cooling function gives: 



A(p,P) Ao.: 



^[e2-ei]n-eip[C2-Cl]n-Cl 



which corresponds to a superposition of an infinity of processes. Introducing this expression into 
once more the integral representation of Gauss hypergeometric function leads to an implicit solution 



(3.3) 

and using 
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(3.4) 



where e„ = ei - [e2 - ei]n, Cn = Ci - [C2 - Ci]n and Ao,„ = ([Ao^2/Aoa]"/Aoa)/5o [Po^^s] ^" We can see this 
solution as the onc-proccss solution with a correction due to the presence of a second process. 



4 Analytical evaluation of fundamental physical quantities 

We will evaluate two fundamental quantities which are the thickness (which corresponds to Xs) of the post-shock 
medium and the accreted column density (that we note S). From the previous results, it is easy to show that: 
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= A[7, A/, e,C] X Vstcooi 



(7 + l)M2(l + l/[7M2])_ 

We can apply this result to magnetic cataclysmic variables in which Vs = \J 2GMuid / Rwd where G, M^d and Rwd 
are respectively the gravitational constant, the white dwarf mass and the white dwarf radius. Furthermore, the 
accretion rate m is given by m = p^Vs- The introduction of these relations in equation (j4.ip with assumptions 
of strong shock {M — *■ 00) and 7 ~ 5/3 provides a simple relation between Xs, m, M^d and Rwd'- 
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(4.2) 



This relation is consistent with Wu et al. (1994) results. Now we can evaluate the accreted column. By 



Table 1. Satured values (limit M — + 00) of the H-parameter for different cooling processes 



(e,C) /7 


5/3 


7/5 


4/3 


1.1 


1.06 


( -2.35, 2.50) 


1.991 


2.894 


3.406 


11.387 


19.022 


( 1.50, 0.50) 


6.967 


10.296 


11.961 


35.288 


57.509 


( 2.50, -0.50) 


5.869 


8.651 


10.047 


29.629 


48.292 


( 2.00, 0.00) 


6.286 


9.273 


10.769 


31.756 


51.752 


( 1.00, 1.00) 


8.302 


12.313 


14.316 


42.327 


68.994 


( 3.00, -1.00) 


5.585 


8.232 


9.560 


28.208 


45.983 



definition, we have S = (/g^"' p{x).dx^ / (Jp^° pQ.dx^. This estimation is achieved for the former five analytical 
solutions by Laming (2004). Our results are presented in figure [Hand numerical results in table [T] These results 
generalize Laming ones for any cooling function. We find the trivial result which is that, for a given radiating 
process, smaller the polytropic index is, larger accreted matter is. 
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Fig. 2. Representation of tiie accreted column density versus Macii number for bremsstralilung cooling for several values 
of polytropic index. We observe the saturation which is comptatible with values of table 1. 



5 Conclusion 

In this paper we present the general analytical solution to the problem of astrophysical radiating shock in the 
case where we have one- or two-cooling processes. These results are very important to validate the input density 
profile which is at the basis of all stability studies. Our solutions are totally consistent with previous evaluation 
of several physical quantities. Moreover it makes possible to predict more general results. This study can be 
considered as the counterpart of studies relative to optically thick radiative shocks (Bouquet et al. 2000). 
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